Gravitomagnetic time-varying effects on the motion of a test particle 
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We study the effects of a time-varying gravitomagnetic field on the motion of test particles. 
Starting from recent results, we consider the gravitomagnetic field of a source whose spin angular 
momentum has a linearly time-varying magnitude. The acceleration due to such a time-varying 
gravitomagnetic field is considered as a perturbation of the Newtonian motion, and we explicitly 
7-H ■ evaluate the effects of this perturbation on the Keplerian elements of a closed orbit. The theoretical 

predictions are compared with actual astronomical and astrophysical scenarios, both in the solar 
system and in binary pulsars systems, in order to evaluate the impact of these effects on real systems. 
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' I. INTRODUCTION 

7— i ; 

i—,' The similarity between Newton's law of gravitation and Coulomb's law of electricity has been largely investigated 
since the nineteenth century, and many authors suggested that, in analogy with the electromagnetic case, the motion 
i 3 = ' of masses could produce a gravitomagnetic (GM) field [l|. Actually, in General Relativity (GR) a GM field is indeed 
generated by the mass current: more in general, any theory that combines Newtonian gravity with Lorentz invariance 
must include a GM field. It is well known that the field equations of GR, in linear post-newtonian approximation, can 
be written in terms of Maxwell equations for the gravitoelectromagnetic (GEM) fields @, [1] (for some limitations of 
£f) . this analogy see [|J , while for a different approach to gravitoelectromagnetism based on the use of tidal tensors in GR, 
' see @, @|). In this context, the role of the gravitational induction has been recently investigated Q, and it has been 
shown that it is possible to describe gravitational induction in analogy with the Faraday-Lenz law of electromagnetism. 
Furthermore, for stationary GEM fields, the spatial components of the geodesic equation reduce, up to first order in 
v/c, to the Lorentz force law, while for general time- varying fields the corresponding expression is more complicated. 
Time-varying GEM fields can be interesting in actual astrophysical events: for instance when a star, or a planet, 
loses its spin angular momentum, it generates a time-varying GM field, whose effects on the neighboring objects can 
be evaluated. As an example, one may consider the decrease of the rotation rate of the Earth (due to tidal effects): 
this corresponds to a loss of angular momentum which, in turn, produces a time-dependent GM field; however, this 
variability is so small that is not expected to influence experiments aimed at the measurement of the GM field of the 
Earth, such as GP-B @,||. 

In this paper we are concerned with the study of the effects of time-variability of the GM field on test particles. 
Starting from the results obtained in Q, we consider a non-stationary spacetime where the magnitude of the spin 
angular momentum of the source of the gravitational field varies linearly with time (the case of the variation of the 
direction of the spin angular momentum for the effects of external torques was considered in II Oil . while the instability 
of closed Keplerian orbits induced by a time- varying gravitomagnetic field was studied in |11|). In particular, we 
focus on the equation of motion of a test particle, where we consider the relativistic linear GM contribution as a 
perturbation of the Newtonian motion. The effects of such a perturbation on the Keplerian elements of a closed orbit 
are evaluated, and then compared with some actual astronomical and astrophysical events, in order to estimate the 
impact of these effects on real systems. 
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II. MOTION IN A NON-STATIONARY GRAVITO-MAGNETIC SPACETIME 

The solutions of the GR field equations, for a localized slowly rotating source, in linear approximation give raise to 
the following spacetime metric (see Q) 

ds 2 = -c 2 (\ - 2^J dt 2 - ^(A • dx)dt + (\ + 2 ^dx'dxJ , (1) 

where the GEM potentials can be expressed as 

GM G S x r 
* = , A = — , (2) 

r c r° 

in terms of the mass M and the spin angular momentum S of the source; in what follows the Cartesian coordinates 
{x,y,z} are choosen in such a way that the angular momentum is parallel to the axis z: S || e z , where S is the 
unit vector in the direction of S; furthermore, we use the notation |r| = r. In order to solve the field equations the 
transverse gauge has been imposed, which, in terms of the GEM potentials reads 

f + v-(i A )=o. 

In particular, in 7] GEM potentials in the form 



T GA*t G r _ ~ , S x r , . 

* = , A = - [S + Sx (t-t ) — 3- , (4) 

r c r° 

have been considered, where the magnitude of the angular momentum varies linearly with time 1 S = S(t — to) = 
[Sq + Si (t — to)}, and it has been shown that the metric (HJ, with these potentials, represents a solution of the 
linearized GR field equations, in the linear GEM approach, in which all terms of 0(c -4 ) in the metric tensor are 
neglected. Of course, this time-varying solution is consistent with the linear GEM approach if the time variation of 
the magnitude of the angular momentum does not break the validity of the linear approximation 2 ; moreover, the 
instants of time at which the temporal variation starts and ends are not considered. 

Actually, it is interesting to point out that the effective source for the spacetime metric ([T]) with the GEM 
potentials in the form ((4]) is localized and divergenceless, therefore its volume integral vanishes, even though stresses 
distributed throughout space arise, which fall off as r~ 3 for r — > oo: however this fact does not affect the viability of 
the linear GEM approach, since the effects of these stresses on the metric tensor are of 0(c -4 ) (see the discussion in 
S and 0). 

The spatial components of the geodesic equation of the metric ((T|) , with the GEM potentials (U]) up to the dominant 
gravitomagnetic terms are (see [3]) 

dw GMr 2GSxr 2 n 

= — — 3 vxB ' ( 5 ) 

c r° c 



where the gravito-magnetic field is 



B= -(S +S lt )l[3(S-r)r-S] (6) 



The geodesic equation ([S]) can be written in the form 



In order to carry out the perturbative analysis of the motion, we have chosen the constant So in such a way that it corresponds to the 
magnitude of the angular moment at t = to, i.e. when the particle is at the periastron of its unperturbed orbit (see below). 
2 In other words, we suppose to work within a time interval At for which 2|A| <C c 2 , so that the gravitomagnetic term (i) exists and (ii) 
can be considered as a perturbation of the fiat spacetime metric. 
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where W can be thought of as an acceleration perturbing the Newtonian motion. We set 

W = fi+f 2 , (8) 

where 

2GSxr 

f i = ^^' (9) 

and 

2 

f 2 = — v x B. (10) 
c 

The latter term is nothing but the gravitomagnetic Lorentz acceleration, while the former is an azinmthal acceleration 
due to the time variability of the angular momentum S. 

III. PERTURBATION OF THE NEWTONIAN MOTION 

In order to calculate the impact of the acceleration (|8|) on the orbit of a test particle, let us project it onto the 
radial (e r ), transverse (e t ) and normal (e„) directions of the co-moving frame picked out by the three unit vectors. 
To this end let us briefly describe the geometry of the problem we are dealing with: we consider a first set of Cartesian 
coordinates {a;, y, z}, whose origin coincides with the source of the gravitational field; in particular, in this coordinates 
we have S || e z . Then, we introduce another set of Cartesian coordinates {X, Y, Z}, with the same origin: the X axis 
is directed along the ascending node, the Z axis is perpendicular to the orbital plane. The angle between the x and 
X axis is fl, the longitude of the ascending node, while the angle between the z and Z axis is /, the inclination of the 
orbital plane. To summarize, we write the relations between the basis vectors 



ex — cos £1 e x + sin f2 e y , 

ey = — cos / sin £1 e x + cos / cos f2 e y + sin / e z , (11) 
ez = sin I sin O e x — sin / cos Q e y + cos / e z . 

Let r = r(f)e r be the Keplerian orbit. The radial unit vector e r in the plane of motion is described by 

e r — cos u ex ~t~ sin u ey , (12) 

where u = oj + /, u) is the argument of periastron and / is the true anomaly. We remember that the unperturbed 
closed orbit, i.e. the orbit under the Newtonian acceleration in eq. (JT]), is an ellipse described by 

r(f) - (13) 
1 + e cos / 

where a is the semi-major axis and e is the eccentricity. 
We may write the velocity in the form 

[(— sinw — esinw) ex + (ecosw + cosu) ey] , (14) 



where n — ^ GM/a 3 is the un-perturbed Keplerian mean motion related to the un-perturbed Keplerian orbital period 
by Ph = 27r/n. Now, on using eqs. (fT2 j) - (fil| we can calculate the expression of the perturbing accelerations . ([TU)) 
along the unperturbed orbit. In particular, we are interested on the components of these accelerations along the 
radial, tangential and normal (to the plane) directions, which, if we choose without loss of generality, O = 0, are 
defined by the following unit vectors 

e r = cos u e x + cos / sin u e y + sin / sin u e z , (15) 
e t = — sin u e x + cos / cos u e y + sin / cos u e z (16) 
e n = — sin/ e y + cos/ e z , (17) 
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In order to take into account the explicit time dependence, it is useful to introduce the eccentric anomaly E, which 
is related to the true anomaly / by the following relations: 



. cosE-e Vl - e 2 sinE 

cosf = - -, sm/=— — . (18) 

1 — e cos E 1 — e cos E 

Furthermore, in terms of the eccentric anomaly we may write the unperturbed ellipse in the form 

r = a(l - ecosE), (19) 

the time dependence can be expressed in terms of the eccentric anomaly thanks to 

E — e sin E 

t-t = , (20) 

n 

where to corresponds to the passage at the periastron. 

That being said, we are now able to explicitly calculate the components of the perturbing acceleration W, in the 
form 

W r = (f 1 +f 2 )-e r , 

W t = (fx + f 2 )-e t , (21) 
W n = (fi + f 2 )-e„, 



and we obtain 



= 2G C q S IVT^ nS(E) 

c 2 a 2 [1-ecosE] 4 
2Gcos/ Si — nS(E)esinE — 25iecos£' + S^e 2 cos 2 E 
c 2 a 2 [1-ecosE 1 ] 



2Gsin/ I sinwsinE (l — e 2 ) + coswVl — e 2 \e — cos -El ] , , 

W n = Si < o } + (24) 

c 2 a 2 VT^ \ [1-ecosE] 3 J ' 



2Gsin/ „ /7 ^n I 2 sin uj cos E + 2 cos u sin E vl — 



iS(E) 



2 



c 2 a 2 \/l — e 2 1 [1 — ecosE] 

2Gsin/ I 3 sin u; — sin w cos 2 E — cos u> sin E cos E V 1 — e 2 I 

nS(E) < = > e + 

c 2 a 2 vT^ 'I [1-ecosE] 5 J 



2Gsin/ I 3 cos w sin E\/l — e 2 + 2 sin u> cos E 



™S(E) 77 775 e + 



c 2 a 2 \/l — e 2 I [1 — ecosE] 

2Gsin/ I sinw cos 2 E — 3 sinw I ^ 

c 2 a 2 ^Y^? ^ '\ [1-ecosE] 5 J 



where 

S(E) = So + ^-^^. (25) 
n n 

The components (|22II24[) of the perturbing acceleration must be inserted into the right-hand-side of the Gauss 
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equations for the variations of the Keplerian orbital elements: 

2 



da 
lit 
de 
~dt 
dl 
~dl 

dn 

~dt 

du) 

~dt 

dM 
~dT 



eW r sin f + W t (-) 



W r sin f + W t 



1 



cos / + - 1 



na 



i 



(-) cos(/ + w), 
n (-) sin(/ + w), 



na sin /yi — e 2 

-VF r cos/ + Ty t ( 1 



P 



sin/ 



COS/' 



dt' 



n W r 

na 



e 



,, , dco T df2 
2 ' — + cos/— 



(26) 
(27) 
(28) 
(29) 

(30) 
(31) 



where A4 is mean anomaly of the orbit of the test particle and p = a(l — e 2 ) is the semi-latus rectum. In order 
to obtain the secular effects induced by the perturbing acceleration we must evaluate the Gauss equations onto the 
unperturbed Keplerian ellipse ([T9| . taking into account eqs. (fl8|) and average 3 them over one orbital period Pb of the 
test particle by means of 



dt 



1 — e cos E 
2^ 



dE. 



(32) 



Actually we are interested in the explicit calculation of the secular variations for orbits that have small eccentricities, 
as a consequence we neglect terms of order 0(e 2 ). What we get is: 



4 cos/, 
GSx (2 + e) 



cos/, 







G5i (1 + e)cos 2 a; 
-2 ^—3 sin / 



= 2- 



G 



— sin 2lj (1 + e) — tt 



M) = 



-2- 
0. 



G 



-S 1 ! — sin 2lu — 37r 
' 2 



SS'on > cos/ 



(33) 
(34) 
(35) 
(36) 

(37) 
(38) 



We point out that the variations of the semimajor axis, the eccentricity and the inclination are pure effects of the 
time variability of the angular momentum, since they are proportional to 5*i . On the other hand, in the variation of 
the longitude of the ascending node (f36|) and the argument of the periastron (|37)) there is a term proportional to So, 
which is due to the usual (i.e. static) gravitomagnetic field. 

From (|36|) and (|37| we can get the secular variation of the longitude of pericentre w = Q cos / + u: 



(tb) 



G 



Si 



-esin2w + 2tt 



+ 2Squ > cos/ 



(39) 



We used / + (!; = / because over one orbital revolution the pericentre lj can be assumed constant. 
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IV. COMPARISON WITH THE OBSERVATIONS 



Here we will use (|33|) - (l38l) to make contact with observations in real systems. Let us start from the system Earth- 
LAGEOS; indeed, the secular increase of the diurnal rotation period P 

— = 3 x 1(T 10 yr" 1 (40) 
induces a linear decrease of the Earth's angular momentum [l2j 

S = aMR 2 0^ = 5.86 x 10 33 kg m 2 s" 1 (41) 

of 

Si = -5.6 x 10 16 kg m 2 s" 2 . (42) 

For the LAGEOS satellite (a = 12270 km, / = 110 deg), (|4"2")l and (|33p yield a decrease of its semimajor axis of only 

(a) = -2 x 10~ 14 m yr" 1 , (43) 

which is completely negligible with respect to the observed secular decrease of the semimajor axis of LAGEOS of —0.4 
m yr" 1 [13]. 

Another evaluation of the impact of the variation of these time-depending gravitomagnetic effects can be done by 
considering the double-pulsar system PSR J0737-3039 ([II El). 

In order to give a rough estimate of the effect, we consider the gravitomagnetic field of PSR J0737-3039B, acting 
as a perturbation of the motion of the companion star PSR J0737-3039A. On setting S = =ffl, where 1 is moment of 
inertia of PSR J0737-3039B and P its period, we may write 

l-i 

and, hence 

S a = -2nl^ (45) 

In order to give a rough estimate, we can take for the moment inertia the value X = 10 38 kg m 2 (see e.g. [H[); on 
considering the measured parameters of the system PSR J0737-3039 [H[ , we obtain 

S x = -7.8 x 10 22 kg m 2 s~ 2 . (46) 

On using (|33p . where we take cos / = 1 to consider the most favorable condition, we get a variation of the semimajor 
axis of 

(a) = 4.7 x 10~ n m yr" 1 , (47) 
which is negligible, since the uncertainty on a is about 8.1 xlO 5 m (see [l7j l 



V. CONCLUSIONS 



We have studied the motion of a test particle in the gravitomagnetic field of a source whose spin angular momentum 
has a time- varying magnitude, starting from a recent solution of the GR field equations for a localized slowly rotating 
source in linear approximation. We have treated the gravitomagnetic acceleration as a perturbation of the Newtonian 
motion, and we have explicitly obtained the secular variations of the Keplerian orbital elements due to this acceleration. 
Celestial bodies generally have variable angular momenta, so we have focused on some real systems, in order to evaluate 
the impact of these time-depending gravitomagnetic effects. In particular, we have considered the LAGEOS satellites 
orbiting the Earth and the double pulsar system PSR J0737-3039. In both cases, we have evaluated the magnitude 
of the secular variation of the semimajor axis, and we have shown that the impact of the time variation of the spin 
angular momentum of the gravitational source is negligibly small, thus suggesting that it can be neglected in actual 
astronomical and astrophysical scenarios. 
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